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Abstract Analysis of the secondary structures of mRNAs which encode mature peptides shows 
that the location of each codon in mRNA secondary structure has a trend, which appears to be in agree- 
ment with the conformational property of the corresponding amino acid to some extent. Most of the 
codons that encode hydrophobic amino acids are located in stable stem regions of mRNA secondary 
structures, and vice versa, most of the codons that encode hydrophilic amino acids are located in flexible 
loop regions. This result supports the recent conclusion that there may be the information transfer be- 
tween the three-dimensional structures of mRNA and the encoded protein. 
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Since triplet codons were determined, a lot of studies have been made on the structure of 
the genetic code and its function. As far as the structure of code is concerned, RNY pattern 
was proposed (where R denoting purine, Y pyrimidine, and N any nucleotide) (Crick, 1968; 
Eigen et al ., 1979) for a comma-less code based on physicochemical properties of nucleic acids 
and optimal codon-anticodon interactions, which was also supported by the statistical analysis in 
present rRNAs (Eigen, 1985). It has been noticed that RNY patterns are abundant in current 
genes. This is often interpreted as the molecular track of a primeval comma-less code. Howev- ] 
er, Wong et af. (1986) argued on the basis of mutational rate consideration that observed 
RNY preferences should not have survived by random mutational events without selective inter- 
ference. 

The results mentioned above only involve the distribution of codons in DNA single strand. 
In fact, codon distribution also shows some characteristics connecting with biological function in 
double strands. Analyzed the codon distributions of sense and antisense strands in a polynu- 
cleotide double strand, Konecny et af. (1993, 1995) argued in quantitative terms that the ge- 
netic code is "adapted" to simultaneous coding of two proteins in both complementary nucleic 
acids, i.e. informational adaptation. For instance, silent mutations in the sense strand corre- 
spond to conservative ones in the antisense. [t is also not difficult to see that codons for hy- 
drophobic amino acids in one strand are complementary to codons for the hydrophilic amino 
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acids in another, and codons for charged (hydrophilic) amino acids are prevalently comple- 
mented by codons for uncharged. This structure of genetic codons results from the functional 
pseudo-symmetry of the triplets (Konecny, 1995). 

mRNA is the direct template of protein synthesis, so it should be more significant to study 
the corresponding relation between mRNA sequence and the encoded protein because of the ex- 
istence of the non-coding intervening nucleotide sequence in eukaryotic DNA. RNA molecule 
can fold back onto itself in the structure stabilized by hydrogen bonds between complementary 
bases. The biological functions of RNA molecules are always associated with their fold struc- 
tures. Such a fold structure, i.e.secondary structure, contains stems {base-pairing or double 
regions} and loops (single-stranded regions). The stem of RNA secondary structure is similar 
to a DNA double strand, except G pairing with U in RNA. Some elementary studies on codon 
distributions in RNA secondary structure have been made (Hasegawa et al., 1979; White et 
al., 1972; Konecny, 1995), which reveal some structural features of codons and the encoded 
amino acids. However, these work only considered the situation of "triplet-codon pairing" ac- 
tually, namely the three bases of a codon in one strand of a stem just pair with the three bases 
of another codon in the complementary strand. So it was concluded that RNY pattern codons 
are mainly situated in stem regions and the silent mutations in one strand of a stem correspond 
to the conservative ones in the complementary strand, similar to those in DNA. In fact, the 
stems of an RNA secondary structure do not always be paired by "triplet-codon". However, 
they are paired by base complementary. Therefore more detailed studies are needed about the 
distribution of codons, the influences of silent mutations and other structural features in RNA 
secondary structure. 

Recently, we have made a comparative study on the secondary structures of some mRNAs 
and the three-dimensional structures of the encoded proteins. Some positive correlation has been 
found (Zhang et a£. , 1997). This result suggests that mRNA perhaps carries three-dimension- 
al genetic information (spatial structural information) besides one-dimensional linear informa- 
tion. In order to further understand how the three-dimensional information is contained in an 
mRNA secondary structure and how it is transferred to protein, we analyzed the secondary 
structures of 86 mRNAs which encode the mature peptides. Some distribution features of the 
triplet codons in mRNA secondary structures were found. They appear to be in agreement with 
the conformational properties of the coded amino acids. Here the conformational property of an 
amino acid means the conformation in protein taken by the amino acid and its backbone. This 
result is in favor of the conclusion of the correlation between the three-dimensional structures of 
mRNA and protein. 


1 mRNA samples and their secondary structures 


The predominantly used thermodynamic folding algorithm for the secondary structures of 
all kinds of RNAs was originally proposed by Zuker et af. (1981). The prediction accuracy of 
Zuker's algorithm is higher as far as the current knowledge, especially to smaller RNA 
molecules. For example, the accuracy for (RNA is 96% , and for 5S rRNA is 88% (Jaeger et 
al., 1989; Zuker, 1989; Zuker et af., 1991). It is obvious that the prediction accuracy is 
higher as RNA is smaller. We also used Zuker’s algorithm and PCFOLD program to predict 
the secondary structure of mRNA. PCFOLD program has a limitation of nucleotide length< 
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425, so we have to choose mRNA samples according to this limit. 

mRNA samples were mainly taken from EMBL nucleotide sequence database released 37 
December 1993, involving mammals, invertebrates, vertebrates, synthetics and NCBI back- 
bone. We only considered the mRNA sequences which encode the mature peptides and whose 
nucleotide lengths are less than 425, so 67 samples have been collected in this database. In ad- 
dition, we have also collected 19 samples in GenBank database released by December 1996, and 
they have been used to make the comparative analysis of the three-dimensional structures with 
their corresponding proteins. The nucleotide lengths of about 70% of the 86 samples are less 
than 300. Therefore, the prediction accuracy of their secondary structures should be higher. 


2 Loop codons and stem codons in an mRNA secondary structure 


Every mRNA sample was folded with PCFOLD program to obtain its secondary structure. 
l All triplet codons of every sample were sorted into stem codons and 


3 e ere, loop codons according to their locations in the mRNA secondary 
X AG structure. By a stem codon we mean that its three bases are all 
en paired and there is no any bulge base in the complementary region. 

aU All other triplet codons are called loop codons. It should be men- 

Ce tioned that the codons whose bases are paired completely but with 

co ^g bulge bases in complementary (such as the codon UCC at 3' end in 
GU. Fig.1) are defined as loop codons, because RNA backbone occurs 


Fig.1 A segmet of the bending at this location besides the normal propeller-like twist. 
mRNA secondary stu- From our study of the correspondence between the three-dimensional 
cture of LEXG structures of mRNA and protein, loop codon regions in mRNA sec- 
ondary structure are generally related with the turns or the connect- 

ing peptides of the protein (Zhang et al., 1997; Liu et a1. , 1998). 


3 The distributions of all kinds of triplet codons in mRNA secondary structures 


Eigen e£ al. (1985) grouped triplets in 5S rRNAs into RNY, RNR, YNY and YNR, 
and found their occurring frequencies as following; RNY > RNR > YNY > YNR. We found 
that this relation is primarily remained in mRNAs after analyzing 86 samples, see Table 1, 
where we grouped codons in more detail than Eigen. But to different species, there are some 
differences. It can be seen from Table 1 that the occurring frequencies of RNR pattern are 
higher than RNY’ s in mammals and vertebrates; the frequency of YNY is higher than RNY’ s 
in vertebrates; and the frequency of YNR>YNY in synthetics. These differences maybe result 
from the evolution of species. 

The distributions of all patterns of the triplets in the nucleotide sequences have supplied a 
lot of meaningful information, which are helpful for us to understand the origin and the struc- 
tures of genes. However, considering the actual state (fold structure) of RNA molecules in liv- 
ing cells, the location distribution of every pattern of the triplet codons in the secondary struc- 
tures of mRNAs should also be considerable. Our analysis shows that some meaningful distribu- 
tion features do exist, which seem to be just related to the rules establishing RNA secondary 
structures. Table 2 lists the ratio of the two numbers for every pattern of codons locating in the 
loop regions and in the stem regions respectively, which is called the loop-stem ratio. The loop- 
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stem ratio, to a certain extent, reflects the probability for a codon locating in the stem region or 
in the loop region. The bigger the loop-stem ratio is, the bigger the probability for a codon lo- 
cating in the loop regions is; or vice versa. According to our above definition of a loop-codon 
and a stem-codon, the loop-stem ratio of all triplet codons in an RNA secondary structure is 
generally bigger than one. To our 86 samples, it is 2.34 on the average. It can be seen from 
Table 2 that the ratios have some differences to different species. The loop-stem ratio of RRR 
pattern is highest in the overall result, next is RRY. The loop-stem ratios of RRR in mammals 
and in synthetic are also higher. However, the loop-stem ratio of RRR is lower than the loop- 
stem ratios of YYY and RYY in invertebrates; YYY's is highest and those of RRY and YRY 
are all lower than the average. 

From the statistical results, the loop-stem ratio of the codon seems to be influenced by its 
occurring frequency to some extent. In overall, the loop-stem ratios of the triplets primarily in- 
crease with the increment of their occurring frequencies (Fig.2). But the increment rate is 
small when the frequencies are less than 15.74. The increment rate of loop-stem ratio becomes 
big when the frequencies > 15.74. That is, the loop-stem ratios of RRY, RYY and RRR are 
associaed with their frequencies closely. The situation in mammals is almost similar to that in 
overall. Most of the curve of vertebrates is in agreement with that of overall. However, in in- 
vertebrates and synthetics, there are bigger deviations compared with other species in the mid- 
dle parts, especially in synthetics (Fig.2). 

4 The locations of triplet codons in mRNA secondary structures and the confor- 
mational properties of amino acids 

Synonymous codon usage biases to various secondary structures in proteins has been detect- 
ed (given in elsewhere). Therefore, it should be needed to consider the effect of codon usage 
in studying protein folding, not only to consider amino acids itself. Here, combined with the 
conformational properties of amino acids, we investigate the distributions of codons in mRNA 
secondary structures, 61 triplet codons are grouped into 5 sets according to the loop-stem ratios 
of individual codons (Table 3). Generally, if the loop-stem ratio of a codon is lower, the codon 
is frequently located in stem region of mRNA secondary structure; inversely, if the loop-stem 


Table 1 Occurring frequencies (%) of all kinds of codons in the mRNAs of different species 
Species RRY RYY RRR RYR YYY YRY YYR YRR 
Over? ^ 16.74 15.74 17.57 10.58 41,93 11.88 9.2 6.35 
MAM 16.17 13,61 22.86 11.28 12.59 8.71 7.78 7.00 

INV 16.79 17.44 13.76 11.6 13.00 14.84 7.58 5.53 
VRT 13.57 13.31 17.29 11.86 12.56 17.78 8.27 5.36 
SYN 16.76 14.86 17.71 8.19 10.09 8.76 14.48 9.14 
(MAM denotes mammals, INV invertebrates. VRT vertebrates and SYN synthetics. 
Table 2 Ratio of the numbers for every pattern of codons locating in loops and In stems (loop-stem ratio) 


Species RRY RYY RER RYR YYY YRY YYR YRR Average 
Overall 2.59 2.19 3.54 2.07 2.09 2.07 2.30 1.89 2.34 














MAM 3.33 1.69 4.25 2.02 1.61 2.61 2.03 2.21 2.47 
INV 1.98 3.60 2.85 3.64 4 1.40 2.68 2.64 2.85 
VRT 2.52 2.64 2.97 1.81 2.16 1.97 2.20 1.36 2.21 


SYN 3.89 1.79 5.64 2.07 3.82 2.29 4.43 2.43 3.30 
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Fig. 2 The relationship of loop-stem ratios of all kinds of pattern codons 
and their occurring frequencies 


ratio of a codon is higher, the codon is frequently located in loop region. 

At first, putting the first two sets together, we analyze the case of L/S<2.34 (the aver- 
age value of all samples, Table 2). It has been found that most of the codons with L/S«:2.34 
encode the amino acids that are hydrophobic or prone to be inside proteins. It is well known 
that Phe, Leu, Val, Ala and Gly are hydrophobic amino acids and that Trp has the hydropho- 
bic trend although it is of polarity (Lai et al., 1993). It can be seen from Table 3 that the 
amount of their codons come to about 61% in this set of codons. The loop-stem ratio of Trp is 
the lowest (0.83), i.e. the frequency is high for its bases occurring in stem regions. In addi- 
tion, Glu generally forms salt bridge and Cys forms disulfide bond inside proteins. Therefore, 
summed up, about 7196 of the codons with L/S«:2.34 encode the amino acids with hydropho- 
bicity or the trend situated inside proteins. Arg, Ser and Gln are hydrophilic amino acids; their 
codons are only about 29% in this set of codons. However, it should be noted that Arg has hy- 
drophobic trend although it is charged (Lai et af., 1993). So the proportion of codons coding 
for the amino acids with hydrophobicity or the trend situated inside proteins amounts to 86%. 

The situation is just reverse in the set of codons with L/S>4, i.e. most of these codons 
encode the amino acids that are hydrophilic or prone to be located at the surfaces of proteins. 
About 85% of this set of codons code for Asn, Thr, Gln, Tyr, His, Lys and Arg, which are 
of polarity or charged (i.e. hydrophilic). The codons of Leu, Ile are only about 15%. 

In the set of codons with 2.34<L/S<3, codons UUA, CUU, AUU, AUC, CCC, CCA 
and CCG code for hydrophobic amino acids Leu, Ile and Pro. The codons coding for hy- 
drophilic amino acids are AGU, AGC (Ser), ACG (Thr), GAC (Asp) and CGA (Arg) 
(Asp and Arg are polar). In the set of codons 3< L/S<4, the codons coding for hydrophobic 
amino acids are CUC (Leu), AUG (Met) and CCU (Pro); and the codons coding for hy- 
drophilic amino acids are UCA (Ser), UAU (Tyr), CAC (His), GAU (Asp) and GAA 
{Glu) (Asp and Giu are polar). It can be easily seen that half of the codons whose loop-stem 
ratios are between 2.34 and 4 code for hydrophobic amino acids, and another half of them code 
for hydrophilic or polar amino acids. Of course, Arg and Glu have the trend situated inside pro- 
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Table 3 Codons and corresponding amino acids classified 
according to the loop-stem ratios of codons 
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ratio (LAS) 
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GGU Gly 3«L/S«4 CUC Leu 

CGG Arg AUG Met 

1< L/5«2.34 UUU, UUC Phe UCA Ser 
UUG, CUG Leu CCU Pro 

GUU, GUC, GUA, GUG Val UAU Tyr 

UCU, UCC Ser CAC His 

GCU, GCC, GCA, GCG Als GAU Asp 

CAG Gin GAA Glu 

GAG Glu L/S»4 AAU, AAC Asn 

UGU, UGC Cys CUA Leu 

AGG, CGU, CGC Arg AUA lle 

GGC, GGA, GGG Gly ACU, ACC, ACA Thr 

2.34« L/S«:3 UUA, CUU Leu CAA Gln 
AUU, AUC lle UAC Tyr 

AGU, AGC Ser CAU His 

CCC, CCA, COG Pro AAA, AAG Lys 

ACG Thr Arg 


5 Discussion 


As far as genetic code is concerned, many rules have been found in it, such as the rules of 
the degeneracy. In addition, there are also some meaningful patterns of base distributions relat- 
ed to the properties of amino acids. For examples, the codons whose middle bases are U encode 
hydrophobic amino acids, and the codons whose middle bases are A encode hydrophilic amino 
acids. Some of these rules have been explained from the stereochemical theory. From DNA or 
mRNA linear sequences, the distributions of codons also show some characteristics, such as the 
3-periodicity (Voss, 1992) in coding sequences and the bias of the codon usage (Lid, 1994). 
All of these rules suggest that the assignment and the distribution of codons in genes are not 
random. In another hand, the non-random distribution of codons is also embodied in the three- 
dimensional structures of nucleic acids molecules. This is noticed in the pseudo-symmetry of the 
function in DNA double strands (Konecny, 1995). Mature mRNA is a direct template of pro- 
tein synthesis. It is generally thought that mRNA only provide the linear information to the 
protein until now. However, mRNA has a specific fold structure. The functions of RNAs are 
associated with the three-dimensional structures. 

In the consideration of mRNA secondary structure, generally, a triplet codon can occur at 
any location in an RNA secondary structure. Considering the conventional base-pairing rule of 
RNA (G=C, A=U and G - U) and the base-pairing stacking interaction which influence the 
free energy of the RNA secondary structure, there is a certain location bias for every codon in 
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the secondary structure. We can roughly calculate the probability of every codon occurring in 
stem region or in loop region. Suppose N1N3N3 (N; denotes any base, i 2 1, 2, 3) isa 
triplet codon, only the base-pairing form is concerned, its probability occurring in stem region 
(three bases are paired) can be calculated with 

P, (N,N2N3) =P (N1) P (N2) P (N3) (1) 
where 

-fi if N;=AorC 

P (N) “9, if N;=Gor U 
However, we should note that the interactions of base-pairs are not the same to different base- 
pairs, so are the stacking interactions of the paired bases. For example, G=C pairing is more 
stable than A= U pairing; A= U pairing is more stable than G — U pairing. [n addition, there 
are also some other factors, such as the occurring frequencies of codons, to influence the loca- 
tion of the codons in mRNA secondary structure. So it is necessary to consider these interac- 
tions and the influence factors in order to reflect the probability more accurately. This can be 
achieved by adding certain weights. Thanks to the thermodynamic data available, the weight 
concerning the interactions and the stacking of base-pairing can be determined easily, But the 
weight concerning the other factors is more complex. These problems will be discussed else- 
where in detail. 

According to (1) and combined with base-pairing stacking energy (Turner et al., 
1987), we can roughly estimate the probability of some codons occurring in stem region or in 
loop region. For instance, the codons containing base “A” should occur in loop regions more 
frequently. And the more a codon contains the base "A", the smaller the probability for it lo- 
cating in the stem region is, or the bigger the probability for it locating in the loop region is. 
Similarly, the more a codon contains the base “G”, the bigger the probability for it locating in 
the stem region is. This is in good agreement with the statistical result (Table 3). The loop- 
stem ratios of AAA and AAG (Lys) are 10.15 and 6.09 respectively; the loop-stem ratios of 
AAU and AAC (Asn) are 5.67 and 7.41 respectively; and that of GUU (Trp) is 0.82, etc. 

In RNA secondary structure, the stem region is comparatively stable, inversely, the loop 
region is flexible and has a bigger contacting area with solvent. It is known that the conforma- 
tional property of Lys is flexible and prone to be on the surfaces of proteins. The conformational 
property of Asn is prone to be in the turn regions of proteins. And Trp has the position trend 
inside proteins. So the location distributions of these codons are accordance with the conforma- 
tional trend of the corresponding amino acids. 

Our result shows that the distributions of the codons in mRNA secondary structures are 
not random. These non-random distribution regulations are associated closely with the confor- 
mational properties of the amino acids. This supports our conclusion that there is correlation be- 
tween the three-dimensional structures of mRNAs and proteins. From the viewpoint of infor- 
mation transfer, it is suggested that the conformational properties of amino acids maybe result 
from the distribution features of the codons in mRNA secondary structure. 

It should also be mentioned that there are two particular codons which encode hydrophobic 
amino acids in the case of loop-stem ratio 4. The codon CUA (Leu) and AUA (Ile) have 
higher loop-stem ratios (15 and 4.4 respectively) and they occur in the lowest frequency in 
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their cognate codons (encoding Leu and [le respectively) in our samples. In effect, they just 
belong to the suppressed codons defined by Lió et af. (1994) according to the terms of strong 
and weak bases. It is known that the synonymous codon usage of an amino acid has bias. This 
bias is reflected in the secondary structures of proteins in particular. This result suggests that 
the codon usage bias may be related to the structural need of the protein. More detailed study 
on this aspect will be given in another paper. 
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构 中 的 位 置 有 一 定 的 倾向 性 ， 这 种 倾向 性 似乎 与 相应 氨基 酸 的 构象 性 质 相 一 致 。 大 多 数 编 
码 朴 水 氨基 酸 的 密码 子 位 于 mRNA 二 级 结构 中 较 稳 定 的 鞋 区 ; 反之 ， 大 多 数 编码 亲 水 氨 
基 酸 的 密码 子 位 于 柔性 的 环 区 。 这 个 结果 支持 了 最 近 得 到 的 关于 mRNA 与 蛋白 质 之 间 存 
在 着 三 维 结构 信息 传递 的 结论 。 

XN “三联 体 密码 子 ， 氨 基 酸 ，mRNA 二 级 结构 
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小 能 猫 (Ailurus fulgens) 产 于 喜 玛 拉 雅 - 横断 山脉 地 区 ， 是 地 史 保 存 下 来 的 第 三 纪 子 遗 动物 。 由 
于 长 期 的 生存 竞争 和 自然 选择 结果 ， 使 其 苦 养 生态 位 由 高 质量 、 竞 争 大 的 肉食 转向 低 质量 、 竞 争 小 但 普 
养 成 分 稳定 的 竹 类 。 其 头骨 、 白 货 和 前 掌 都 已 特 化 ， 适 于 咀 咽 和 抓 食 竹 子 。 但 由 于 其 消化 道 仍 是 典型 的 
食肉 型 ， 只 能 从 人 竹子 细胞 的 可 洲 性 物质 中 获取 能 量 。 而 竹子 的 植物 纤维 含量 离 ， 细 胞 可 洲 性 物质 含量 
低 ， 小 能 猜 必 然 要 采取 某 种 优化 的 生态 对 策 来 保证 自身 生存 、 繁 衍 的 能 量 需 要 。 本 研究 旨 在 从 生境 和 食 
物 的 选择 方面 对 该 地 区 小 熊猫 的 生态 对 策 进行 探讨 。 


1 研究 地 点 


高 黎贡山 自然 保护 区 位 于 云南 西部 98"34' 一 98'50' 下 ，24'56' 一 26"09'N， 北 接 青 藏 高 原 东 系 ， 南 北 走 
A, PHRMA RE, RMR SAR 10H LH AMS, OAFHBKE4 AWS, WK 
剧烈 ， 海 找 1090~3916 m. MRA TAME RRMA, RRP AR, SERIE EE, UST, XX 
UA SRA, X 

限于 人 力 和 经 费 ， 经 过 预 研 ， 研 究 区 域 选 在 片 马 了 口 附 近 ， 海 拔 2 800 一 3 400 m 的 30 krè HER 
内 。 该 区 域 地 势 复 杂 陡 峭 ， 气 娱 随 海 找 、 坡 向 变化 较 大 。 站 口 年 均 气 温 TC, € 12.5C (8A), EM 
-0.SC (2 月 )，12 月 至 次 年 3 月 为 积 直 期 。 植被 主要 有 : (1) 温 凉 性 针 叶 林 ,分 布 于 海 找 2 800 一 3 200 
mo 乔木 主 要 是 云南 铁 杉 【Tsuga dumisa). WH (Pinus armandi), RZEK ROM, WARS, + 
木 及 冷杉 ， 成 层 现象 明显 。 林 下 灌 从 有 第 人 竹 、 杜 鹃 ， 主 度 随 林 形 和 地 势 而 蜡 ， 一 般 达 50%, 高 者 可 达 
9096. (ZIRMA, EMR 3 200 一 3 400 m, ELE TT, TENCRUECBSOE XE, 3E ER PUE UD E 
(Abies delavayi )， 下 层 草 本 多 为 县 类 与 禾 本 科 植 物 。 

BESE IR DAR RR TIEN: UEM ir (Fargesia arbiculata )， 分 布 在 海 找 3000 一 3400 m, A 
1.5—2.5 m, WE 5—66 Roe, HRM REATIOM. QARMT CF. papyrifera )， 分 布 于 海拔 1800 
—3200 m, K2—5 m, 密度 4 一 12 k/m. 


2 研究 方法 


沿 山 背 在 东西 坡 海 拔 2 800— 3 400 m 间 分 别 设 4 冯 样 带 , 共 8 条 。 按 一 定时 间 间 陋 ( 约 2 个 月 ) i 
样 带 搜 侯 并 记录 各 样 带 上 小 熊猫 活动 痕迹 〈 采 食 残 迹 、 凌 便 ) 及 相应 生境 情况 (海拔 、 坡 度 、 人 竹子 密度 、 
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